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ABSTRACT: In this paper, a discrete time inventory system with demands is generated by a finite 

number of homogeneous sources. The inter demand times of the primary arrival follows geometric 
distribution. The inventory is replenished according to an (s, S) policy and the lead times are assumed to 
follow a geometric distribution. The demands that occur during stock out period is permitted enter into the 
orbit. These orbiting demands retry for their demand after a random time, which is distributed as 
geometric. These orbiting demands send out signal to get their requirement and the time inten’al between 
two consecutive requests is distributed as geometric and the retrial policy is the multiplicative retrial 
policy. The joint probability distribution of the number of demands in the orbit and the inventory level is 
obtained in steady state case and some important system performance measures are derived in the steady 
state case. 

I. INTRODUCTION 

Inventory models have been considered under continuous review as well as periodic review. In the 
recent past discrete time models have started receiving attention of researchers in the areas of queueing and 
telecommunications [2]. In discrete time setting, it is assumed that the time axis is calibrated into epochs by 
small units and that all the events are deemed to occur only at these epochs. With the advent of fast computing 
devices and efficient transaction reporting facilities, such epochs with small gaps can be conveniently assumed 
so that events can occur at these epochs. 

In the case of inventory modelling under discrete times, the first paper was by Bar -Lev and Perry 
[5], who assumed that demands are non-negative integer value random variables and items have constant life 
times. Lian and Liu [8] developed a discrete time inventory model with geometrically distributed inter -demand 
times, bulk demands and constant life time for items. They assumed (0, S) ordering policy, with instantaneous 
supply which clears any backlog and restores the stock to the maximum capacity S. This assumption helped 
them to have fixed life time for all items. They derived the limiting distribution of inventory level through 
matrix-analytic method. D. L. Minh[10] developed the model of discrete time queue on the finite population 
with inter arrival time follows binomial distribution and general service times with first come first serve basis. 
Abboud [1] studied a discrete inventory model for production inventory systems with machine breakdowns. 
They assumed that the demand and production rates were constant and that the failure and repair times of each 
item were independently distributed as geometric. Lian et al. [9] developed a discrete time inventory system 
with discrete PH-renewal process for (batch) demand time points and assumed discrete- PH-distribution for life 
time of items. They also assumed zero lead time and that unsatisfied demand were completely backlogged. 

I.Atentia and P. Moreno [4] developed a discrete time retrial queueing system with multiplicative 
repeated attempts. The authors used the retrial policy is the multiplicative one which analogous to the linear 
retrial policy in continuous-time[3]. The rest of the paper is organized as follows. In Section 2, we describe the 
problem under consideration and in the next section we model the problem mathematically. The steady-state 
analysis of the model is presented in Section 4 and some important system performance measures are derived in 
Section 5. 

Notation: 

• [A]ij : 

• 0 

• I 

• e 



entry at (i, j) th position of A 
zero vector 
identity matrix 

a column vector of l’s with appropriate dimension 
Kronecker delta function 
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II. PROBLEM FORMULATION 



We consider a discrete time retrial inventory system where the time axis is divided into intervals of 



equal length, called slots (epochs). It is assumed that all system activities (arrivals, retrials and replenishment) 
occur at the slot boundaries, and therefore they may occur at the same time. The maximum capacity of the 
inventory is S. 

• The arrival of demands generated from a homogeneous population of finite size(M). Each demands 
generates the arrival according to a Bernoulli stream with rate 1— p, thus 1— p is the probability that a 
demand arrives at a slot and p is the probability that an arrival does not take place in a slot. When the on- 
hand inventory level is more than one then the arriving demand is satisfied immediately. 

• We adopt (s, S) ordering policy, that is, when the number of available items reaches the value s, an order for 
Q(= S~s > s) items is placed which is delivered after a lead time of geometric distribution with parameter q 
> 0. The condition S - s > s is assumed so that when the supply of an order is received during the stock out 
period, the inventory level would be brought above the reorder level. This provides an opportunity to place 
a reorder when a demand bring the inventory level down to the reorder level. Otherwise the inventory will 
have perpetual stock out. 

• If the arrival finds the inventory level zero, he enters into the orbit of pending demands and he retry after a 
random amount of time. If more than one repeated demand retries at the same slot, any of them is randomly 
selected and the others must go back to the retrial group. The time between two successive repeated 
attempts is geometrically distributed with probability 1 - r 

1 — S " ,k , p ( ((l ] , r 2 )£■[(), 1], l ] , r 2 F \)) given that there are k demands in the orbit. 

Unlike continuous review inventory systems, multiple events such as demand, supply and retrial from the orbit 
may occur between epochs n and n+1, n = 0, 1, 2, . . . . Hence we adopt the following convention: 

If the events such as demand for an item, retrial from the orbit and supply of an order takes place at n (n = 1,2, 
3. . . .), it is assumed that first supply is received then demand occurs and finally retrial from the orbit takes 
place. 



Let Xn denote the number of demands waiting at the orbit, and Ln denote the inventory level at time n. 
From the assumptions made on the input and output processes, it can be shown that the stochastic process (X,L) 
= {(Xn,Ln), n 6 N} is a Discrete Time Markov Chain with state space given by 

E = {(i, j) : i = 0, 1, . . . ,M, j = 0, 1, . . . , S}. 

The transition probability function is defined as for (i, j), (k, 1) £ E, 

p((i, j), (k, 1)) = Pr[Xn+ 1 = k,Ln+l = HXn = i,Ln = j]. 

The transition probability matrix P of this process, 

P = (( p((i,j), (k, 1)) )), (i,j), (k, 1) £ E 

Then the transition probability matrix P can be viewed as. 



III. MODEL DESCRIPTION 



0 1 2 



3 



M -1 M 



0 A 0 B 0 

1 C, A, B . 



P= 2 
3 



C 2 A 2 b 2 



M -1 
M 




where. 
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= 0,1,2, 5 








q(X-p M ) 


j = Q + i- 1 


i = 


= 0,1,2, 5 




[4>]« 


— , 


qp M 


j = i 


i = 


0,1,2, s 






P M 


j = i 


i = 


0,1,2, ,s 








1 -p M 


j = i~ 1 


i = 


s + 1, s + 2, 


5 






0 


otherwise 








For k = 


0, 


l, . . . ,M — 1, 










[B k ],, 




q(\-p M - k ),j = 


0 , i = 0 












0, otherwise 










For k= 


i,2,: 


3 M 














p M ~ k qrA 


j = Q + i 




i = 0,1,2 


5 






1 

§ 

i 

*- 
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j = Q + i 


-1 


i = 0,1,2,.... 


5 






(1 -8 kM )p M - k q 


+d k , M q J = 0 




i = 0, 








(1 — p M k )q 


j = i 




i = 1 




[AJ, 
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p M ~ k qr A 


j = i~ 1 




i = U, 


..5 






(\-p M - k )qrA 


j = * ~ 1 




i = 2,3, 


..5 






M-k k 

P h r 2 


j = i 




i = s + 1, s + 2,. 


s 






(1 - P M ~ k )rA 


j = i ~ 1 
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= s + 1, , 


...s 
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otherwise 












p M ~ k q(\-rA) 


j = 


Q + i 


-1 i = 


0,1,2, 5 






a-p M ~ k ) q a- 


r i r i) J = 


Q + i 


-2 i = 


0,1,2, 






p Mk q{\-rA) 


j = i 


-1 


i — 


1,2,3, , 


[AJ, 


= - 


( \-p M - k )q{\ - 


r Pi) j = * 


-2 


i = 2,3/ 


1 5 






P (1-n r 2 ) 


j = i~ 


1 


i - 5 + 1, 5 + 2, S 






(1 - p M k )(1 _ r l 


4 ) j=i-l 




i - s + 1, 5 + 2 


5 



0 otherwise 



Here we notice that ,all the matrices A k , B k and C k are square matrix of order S+ 1 . 

IV. Calculating Limiting Probabilities 

It can be seen from the structure of P.the homogeneous Markov chain {( X n , L n ), n eNjon the 
finite state space is irreducible. Hence the limiting probability distribution 

x (lJ) = P r [X n = i, L n = j] I X 0 =k,L 0 =l], 

Where Tr' 1 ' 1 ' is the steady state probability for the state (i.j) exists and is independent of the initial state 
(k,l).Letn be the steady state limiting probability vector of P. That is II satisfies 

rip = n 

and 
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Ue = 1 . 

The vector n can be represented by 



(n <0> ,n <l> ,n <2> , 


n <M> ) 




= (n (i ' 0) ,n (U) ,n (i ’ 2) , 


n (,: ’ S) ), for i=0, 1.2, 


M. 



Now the structure of P shows, the model under study is a finite birth death model in the Markovian 
environment. Hence we use the algorithm discussed by Gaver et al.[6] for computing the limiting probability 
vector. For the sake of completeness we provide the algorithm here. 

Algorithm: 

1. Determine recursively the matrix Dn, 0 □ n □ M by using 

Do = A 0 

D; = A, + Q(I - Di-O-lBo, i = 1, 2, . . . ,M. 

2. Solve the system 

n <M> (i -d m ) = o 

3. Compute recursively the vector IT ' ’, / = M —\, 

n <!> = (rr +1) c. +1 (/ - d, y\i = M- 1 ,... 

4. Normalize the vector n , by using 

Ue = \. 



...0 using 
.0 



V. System Performance Measures 

In this section, we derive some system performance measures in the steady-state case. 

5.1 Expected inventory level 

Let lL denote the expected inventory level in the steady-state. Then _i is given by 

M S 

£ = S2>‘' J> 

1=0 j = 1 

5.2 Expected reorder rate 

Let ff denote the expected reorder level in the steady-state. Then _r is given by 



c,.=(l -P M )n 



r (0,s+l) 



+ (1 - P M ~‘ ) + (1 - r x r[ )P M ~‘ ));r (, ' s+1) + (1 - P M { )(1 ■ 



i)7T (U+2) ] 



5.3 Expected number of customer in the orbit 

Let C, {) denote the expected number of customer in the orbit in the steady-state. Then _o is given by 

M MS 

Co = Z«n < ' > e = XX /;r< '’ 7) 

,=1 i = 1 7=0 

5.4 Fraction of successful rate of retrials 

The fraction of successful rate of retrials is given by 

MS M 

^(1 - r^ )7t (l,i) + (1 - 

The successful rate of retried ~ i~\ 

The over all rate of retrial vv 1 /, an 

i-l 7-0 

5.5 Cost function 

The long-run expected cost rate for this model is defined to be 
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C(S,s,M) = c h C i +c s C r +c w C o 

where, 

cs : Setup cost per order 

ch : The inventory carrying cost per unit item per unit time 
cw : Waiting cost of a customer in the orbit per unit time 
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